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THEORY OF WING-BODY DRAG AT SUPERSONIC SPEEDS 
Ey Robert- T. Jones 


At subsonic speeds the pressure drag arising from the thickness of 
the body or wings is negligible so long as the shapes are sufficiently 
well streamlined to avoid flow separation. In that range there exists 
no possibility of either favorable or adverse interference on the pres- 
sure distributions themselves. If one body is so placed as to receive 
a drag from the pressure field of another then the second body is sure 
to receive a corresponding increment of thrust from the first. 

At supersonic speeds this tolerance, which was permitted the 
designer, disappears^ and the drag becomes sensitive to the shape and 
arrangement of the bodies. To be sure, the primEcry factor here is the 
thickness ratio, but nevertheless there exist arrangements in which a 
leirge cancellation of drag occurs. Examples of the latter are; the 
sweptback wlag and the Busemann biplane. 

Recently R, T. Whitcomb (ref. l) has shown how the drag at tran- 
sonic speeds may be reduced to a surprising extent by simply cutting 
out a portion of the fuselage to compensate for the area blocked by 
thf wing. The purpose of the present paper is to discuss some of the 
theoretical aspects of this method of drag reduction and to show how 
the basic idea may be extended to higher speeds in the supersonic range. 

Whitcomb’s deduction of the "area rule" was based on considerations 
of stream tube area and the phenomenon of "choking" — which follow from., 
one-dimens lonal-f low theory. Each individual stream tube of a three- 
dimensional-flow field must obey the law of one-dimensional flow. While 
we cannot actually determine the three-dimensional field on this basis 
alone, nevertheless it provides a good starting point for our thinking. 

The results demonstrate again the effectiveness of basic and simple 
considerations . 

While one-dlmensional-flow theory thus provides a clue to the area 
rule, the necessary principle appears more specifically in the three- 
dimensional-flow theory. Thus, the formulas for wave drag given by 
linear theory, if followed toweird the limit as M approaches 1.0 (from 
above), show that the wave drag of a system of wings and bodies depends 
solely on the longitudinal area distribution of the system as a whole. 

This was first noted by W. D. Hayes in his 19^ thesis (ref. 2). However, 
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because of the limitations of the theory at transonic speeds, this result 
was not thought to be of practical significance. Later G. N. Ward 
(ref. 3), E. W. Graham (ref. U), and others, restricting themselves to 
very narrow shapes, expressed the wave drag in terms of the longitudinal 
area distribution for Mach numbers above 1.0, where the linear theory 
has a better justification. 

It should be noted, however, that both of the problems cited are 
limiting cases of the more general problem of supersonic drag and it 
should be borne in mind that only in certain cases has it been possible 
to reduce the general theoretical formulas to the form of an area rule. 

It can be shown that the flow field about any system of bodies may be 
created by a certain distribution of sources and sinks over the surfaces 
of the bodies. Hayes' formula and the formulas given in reference 5 
relate the drag of such a system to the distribution of these singulari- 
ties. To obtain a formula for the wave drag in terms of area distri- 
butions we have to adopt a simplified relation between the source 
strength and the geometry of the bodies, namely, that the source strength 
is proportional to the normal component of the stream velocity at the 
body surface. There are examples (e.g., Busemann biplanes and ducted 
bodies) for which this assumption is not valid. If, on the other hand, 
we limit ourselves to thin symmetrical wings mounted on vertically sym- 
metrical fuselages, there are indications that a good estimate of the 
wave drag at supersonic speeds can be obtained on the basis of the simpli 
fied relation assumed. 

Following Hayes* method of calculation, we find that at M = 1.0 the 
expression for the wave drag of a system of wings and bodies reduces to 
Karman's well-known formula for the wave drag of a slender body of revo- 
lution,, that is, 

Qxra p+l/2. p + l/Z 

% _>i= — / / S"(x) S"(xj^) log Ix-X 3 _] dxdxj_ 

^-Z/2 

Here S(x) represents the total cross-sectional area intercepted by a 
plane perpendicular to the stream at the station x (see fig.) and 
S"(x) is the second derivative of S with respect to x. Following 
Sears (ref. 6) we may expand S'(x) In a Fourier series and obtain in 
this way a formula for the drag which is completely analogous to the well 
known formula for the induced drag of a wing in terms of its spanwise 
load distribution. Thus, if we write 

X = Z/2 cos cp 


and 


S'(x) = ZAn sin n cp 
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Area distribution given by intersections of Mach planes. 


field it is necessary to superimpose disturbances at all of these angles 
and, in computing the drag, to consider the drag of all the equivalent 
bodies of revolution. The final value of the drag is simply the average 
of the values obtained through a complete rotation of the Mach planes. 

In order to make these statements more specific, we may write the 
equation of one such Mach plane as follows: 


X = x-y' cos 0 - z' sin 0 

where y' = y, z* = v^rF-1 z, and 0 is the angle of rotation 

of the Mach plane. By assigning different values to X while keep- 
ing 0 constant, we obtain a series of parallel planes at the same 
angle 0 around the x axis. By assigning different values to 0 
while keeping X a constant, we obtain a set of planes enveloping that 
Mach cone whose apex lies at the point X = x. 

Selecting a value of 0, we cut through the wing-body system with 
a series of planes corresponding to different values of X. The total 
intercepted area in each plane is then equated to the area intercepted 
by this plane passing through the equivalent body of revolution. If we 
denote the area intercepted obliquely by s(X,0), then the area S(X,0) 
is defined by 
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where ^ is the Mach angle (i.e., sin ji - 1 /m). Thus, S is the area 
intercepted by norma],, planes passing through the equivalent, body of rev^ 
olution on the assumption that this body is slender. Again, we write 


S*(X,e) = ^ S(X,0) = SAn sin nq) 


with 


cos m = — 

Xo 

Here, however, both the length 2 Xq and the shape of the equivalent 
body vary with the angle 0. The drag of each equivalent body of 
revolution, which we may denote by D’(0) is then determined by apply- 
ing Sears’ formula: 


D’(0) 


8 


EnAn^ 


The total drag of the wing-body system is the average of all these values 
between 0 = 0 and 0 = 27t, that is, 

1 

D = — / D'(0) d0 

In general, the coefficients A^ will be functions of the angle 
of projection 0. However, calculation shows that the first two 
coefficients A^^ and Ag are again related in a simple way to the base 
area and the volume v. Thus, 


2 S(Xq) 
It Xq 


Ag - 


k V 

« X^ 


None of the higher coefficients contribute to the base area or volume, 
but they invariably contribute to the drag. 

The rules for obtaining a low wave drag now reduce to the rule that 
each of the equivalent bodies obtained by the oblique projections should 
be as smooth and slender as possible, the "smoothness" again being 
related to an absence of higher harmonics in the series expression 





NACA RM A53Hl8a 


Cl 


7 


« 


\ 


will be sufficient to esti- 
mate roughly the width of 
the transonic zone by con- 
siderations such as those 
given in reference 8. In the 
present case it will be noted 
that agreement with the linear 
theory is reached at Mach num- 
bers above about 1 . 08 , and the 
linear theory clearly shows 
the effect of the modification. 



For further theoretical 
studies of wing-body drag, 
shapes have been selected which 
are especially simple analyti- 
cally, namely, the Sears -Haack 
body and biconvex wings of 
elliptic plan form, having aspect ratios of 2.54 and 0 . 635 . The follow- 
ing figure shows the effect of wing proportions on the variation of wave 
drag with Mach number, both with and without the Whitcomb modification: 


Comparison of Whitcomb’s experiments 
with theory. 






Effect of Whitcomb modification on calculated wave drag. 

In each case the modification has the effect of reducing the wave drag 
to that of the body alone at M = 1.0. In the case of the low-aspect- 
ratio wing this drag reduction remains effective over a considerable 
range of higher Mach numbers. With the higher aspect ratio, however, 
the drag increases sharply at higher ' speeds, so that at M = 1.6 the 
modification nearly doubles the wave drag. 


The rapid increase of drag in the case of the high-aspect-ratio 
wing is, of course, the result of the relatively abrupt curvatures 
introduced into the fuselage lines by the cutout. Such abrupt cutouts 
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are necessarily associated with wings having small fore and aft dimen- 
sions, that is, unswept wings of high aspect ratio. 

These considerations led to the problem of determining a fuselage 
shape for such wings that is better adapted to the higher Mach numbers. 
The first step in this direction is, obviously, simply to lengthen the 
region of the cutout — thus avoiding the rapid increase of drag with 
Mach number. The problem of actually determining the best shape for the 
fuselage cutout at any specified Mach number has been undertaken by 
Harvard Lomax and Max. A. Heaslet at Ames Laboratory. Their solution 
of this problem provides a definite method for determining the distri- 
bution of sources and sinks along the fuselage axis that will achieve 
a minimum value of the drag for a given wing shape at any specified Mach 
mmber. Furthermore, by admitting singularities of higher order — 
quadrupoles, etc., which would distort the rotational symmetry of the 
fuselage, they have been able to show that the wave drag of a wing-body 
system can be reduced, in principle at least, to a minimum value associ- 
ated with the given over-all length and volume of the system, that is, 
to the value for a simple Sears-Haack body containing the whole volume 
of the system.^ 

By adopting our simplified relation between the source strength and 
the body shape, we may describe the result of this theory by a rela- 
tively simple concept, which is illustrated by the figure below. For 



Design of fuselage modification for specified Mach number. 

modifications of the first type, the problem is to determine the 
area to be removed from the fuselage to best compensate for a 

given wing. (See above fig.) Selecting a station along the fuselage 
iThis value is, of course, not an absolute minimum for a given volume 
since, as shown by Ferrari, the wave drag of a body can be reduced to 
zero by special volume distributions (see ref. 9)* 
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axis and a Mach plane passing through this station, we revolve this 
plane around the axis, measuring at each angle 9 the nonnal projec- 
tion, or frontal projection, of the area intercepted where the plane 
cuts through the wing. After plotting these areas against 6 and 
integrating between 0 and 2it, we obtain -,ASf> as the average of the 
values of S^, At any Mach number the total voliime to be subtracted 
from the fuselage is equal to the wing volume. At higher Mach niimbers, 
since the modification extends over a greater length, the area sub- 
tracted at individual cross sections becomes less. 


The figure opposite shows 
the calculated result of design- 
ing the fuselage cutout for a 
specific Mach number, 1.2 in this 
case. The lower ciarve is an 
envelope showing the minimum 
values that can be achieved by 
such a radially symmetric cutout. 



ORIGINAL 


MODIFIED FOR M' 


-MODIFIED FOR M> 1.2 


Effect of modification designed for 
a specified Mach number. 


The figure below shows the \/ / 

magnitude of the gain that is V 

possible by higher order modi- qi _ 

fications of the fuselage shape. ’ 

There are three lower bounds r ENVELOPE 

here, and the symbols Uq, aa, O'- Lr J- J- -u ^ 

etc., attached to them refer to ' " M ' ' ' 

a representation of the fuselage 

shape by singularities of increas- Effect of modification designed fo 
Ingly higher order. The curve a specified Mach number, 

labeled ao is that given on the 

previous figure and shows the maximum effect of radially symmetric modi 
fications. While the fuselage shapes for the other curves have not 
actually been determined, the 
curve labeled Bq + may .04 - 

be thought of as referring _ / ] - j ^ 

to a cutout with an addi- r“'“T 

tional elliptic modification. 03- original 

It will be Interesting to pur- WAVE _ / 

sue this investigation fin-ther \ ENVELOPES 

and ascertain just how the , \ y 

fuselage must be distorted to \ /^o +°2 

cancel the wave drag of the r>i 

wing completely, as indi- ~ / / 

cated by the lowest envelope — — 

curve. Of course, it will be 

necessarv tn .Qt.nT*+ vith a cer- ' 

tain minimum diameter in order 

to preserve a real shape. Envelopes for drag at design Mach 

number. 
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AMES 2x2' WIND TUNNEL 



Experiments on bodies with elliptic wings. 



Calculated values for bodies with elliptic 
wings . 


In order to test 
this theory of deter- 
mining optimum body 
shapes we have started 
a program^ using models 
similar to those inves- 
tigated theoretically. 
Several of these models 
have already been tested 
in the Ames 2- by 2-foot 
wind tunnel, with results 
that agree quite well 
with calculations made 
on the assumptions given 
earlier. Shown here are 
the experimental and 
theoretical curves . The 
aspect ratio of the wing 
in these preliminary 
cases is not sufficiently 
high (AR = 2) to enable 
really striking gains to 
be shown. However, it 
is evident that the 
calculated differences 
are all reproduced in 
the experimental values. 
The experimental series 
will include models 
having higher aspect 
ratios, and we expect 
more significant gains 
to appear. 


There are, of course, examples of wing-body systems which would 
hardly benefit by any change in sha^>e of the fuselage. It is easy to 
decide whether a gain is possible, or worthwhile, by comparing the actual 
wave drag of the system with that of a Sears-Haack body containing the 
over-all volume of the system. In the case of the 63° wing-body combi- 
nation, which has been described in several previous reports, this com- 
parison yields 0.00^1-5 as a lower bound for the wave-drag coefficient 
and 0.005 foT the actual value. In such cases, for which the wave drag 
is initially very low, further reduction by reshaping the fuselage is 
not worthwhile. 


It is clear from the foregoing, however, that appreciable savings 
in drag can be made in many cases by a suitable shaping of the fuselage. 
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Unswept wings of high aspect ratio are benefited most and require the 
most careful consideration of the fuselage shape. 

These new developments illustrate, again, the fact that the dis- 
turbance fields at transonic and supersonic speeds are essentially three- 
dimensional phenomena. It was not long ago that our ideas concerning 
the wing section — which had their origin in the older incompressible 
flow theory — had to be relinquished because of the predominating 
effects of the wing plan form. Now we must learn how to design the wing 
and the fuselage together. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., July 8, 1953 
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